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CALIBRATION TECHNIQUE FOR COINCIDENCE IMAGING SYSTEMS 

Background of the Invention 

The present invention relates to medical imaging 
systems. It finds particular application in the 

calibration of positron emission tomography (PET) or 
coincidence -capable gamma camera medical imaging systems, 



1^ and will be described with particular reference thereto, 

ly It is to be appreciated that the invention may also find 

*!3 application in other types of nuclear cameras as well as 

f Sl other types of diagnostic imaging with two dimensional 

i! 10 rotatable detectors. 

:ff A number of types of radiological events produce 

py two emission photons which travel outward in exactly 

^ opposite directions. For example, a positron-electron 

annihilation event produces such an oppositely directed 
15 photon pair. Coincidence imaging systems take advantage 
of this geometric property in spatially localizing the 
radiological event to a line of response (LOR) . The LOR 
is defined as the line connecting two simultaneous 
radiation detection events, which are presumed to 
20 correspond to detection of the two oppositely directed 
photons. Two detection events are typically judged to be 
simultaneous, or coincident, if both detections occur 
within a preselected coincidence time window. The 
radiological event, e.g. the positron-electron 
25 annihilation, will have occurred at some point along the 
LOR. Ideally, for an approximate point source which 
generates many positron-electron annihilation events 
within a very confined space, the LOR corresponding to 
each annihilation event will pass through the point 
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source, and the point source may thus be located in two- 
dimensional or three-dimensional space, dependent upon the 
geometry of the detector system. Similarly, an extended 
radiological source, such as an organ which has absorbed 
an appropriate administered radiopharmeceutical , may be 
studied by reconstructing the LOR's corresponding to 
electron-positron annihilation events generated by the 
radiopharmeceutical into an image. 

Non-idealities in real detection systems, such 
as mechanical misalignments, degrade image resolution. 
For example, misalignment of the radiation detectors blurs 
the image. The precision with which the radiation 
detectors can be mechanically aligned is limited by the 
weight and bulkiness of the shielded detectors, and the 
alignment precision required for adequate image resolution 
is often not practically attainable by purely mechanical 
procedures. Alignment-related resolution problems are 
particularly acute for positron coincidence imaging using 
a multi-head single photon emission computed tomography 
(SPECT) gamma camera with coincidence circuitry, often 
referred to as a gamma-PET imaging system. The gamma-PET 
provides a versatile and less expensive alternative to a 
dedicated PET imaging system, but the typically large 
heads and rotating gantry complicate mechanical alignment. 

Improved resolution can be obtained by following 
the mechanical alignment with a calibration step. The 
calibration step preferably determines correction factors 
for the positional coordinates of each detector. The 
correction factors may subsequently be applied during data 
acquisition or analysis to improve image resolution. 
Gamma-PET systems have been calibrated using a SPECT data 
acquisition mode, the results of which are applied when 
operating the system in a PET or SPECT imaging mode. This 
approach corrects for the tangential detector head 
coordinate. However, SPECT is not strongly sensitive to 
the radial and orientational detector coordinate 
parameters. An additional disadvantage of calibration in 
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SPECT mode is that subsequent to calibration the detector 
heads are re-positioned into the PET configuration, which 
may introduce additional misalignment not accounted for by 
the prior SPECT calibration. 
5 The present invention contemplates a calibration 

procedure for coincidence imaging systems such as 
gamma-PET and dedicated PET systems, which overcomes the 
above shortcomings and others. 

Summary of the Invention 

10 In accordance with one aspect of the present 

invention, a method for calibrating a coincidence imaging 
system which includes a plurality of radiation detectors 
is disclosed. A plurality of coincidence radiation events 
associated with a point radiation source are measured. 

15 Initial values are assigned for a set of fitting 
parameters. A minimization algorithm is applied, which 
includes calculating lines of response (LOR) based upon 
the fitting parameters and the measured radiation events, 
generating a figure of merit characterizing the apparent 

2 0 size of the point radiation source based upon the LOR's, 
and optimizing the fitting parameters to produce a 
minimized figure of merit. After the minimization, a 
correction factor relating to a positional coordinate of 
a detector is extracted from the fitting parameters. 

2 5 In accordance with another aspect of the present 

invention, a method for imaging using a plurality of 
radiation detectors is disclosed. A plurality of 
coincidence radiation events associated with a point 
radiation source are measured. Initial values are 

3 0 assigned for at least one fitting parameter. Lines of 

response (LOR) are calculated based upon the at least one 
fitting parameter and the measured radiation events. A 
figure of merit is generated that characterizes the 
apparent size of the point radiation source based upon the 
35 LOR's. The at least one fitting parameter is optimized 
using a minimization algorithm which includes iteratively 
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repeating the calculating and generating steps to produce 
a minimized figure of merit. At least one correction 
factor is extracted from the at least one optimized 
fitting parameter. A set of radiation data is acquired 
from an associated subject. The radiation data is 
corrected for camera misalignment by correcting the 
spatial coordinates of the detected radiation events using 
the at least one correction factor. An image 

representation is reconstructed from the corrected 
radiation data. 

Preferably, the at least one fitting parameter 
includes a parameter related to the radial positional 
coordinate of a detector, a parameter related to the 
tangential positional coordinate of a detector, and a 
parameter related to the orientational positional 
coordinate of a detector. For a multiple-head imaging 
system, the fitting parameters preferably include: Ar^ 
i=l to N, where Ar. is a correction for the radial 
coordinate of the ith detector; Atj, j=l to N, where Atj is 
a correction for the tangential coordinate of the jth 
detector; and A9 k , k=2 to N, where A0 k is a correction for 
the orientational coordinate of the )cth detector. 

The figure of merit is preferably generated by 
summing the distance or the square of the distance of 
closest approach of each LOR to a spatial point, in which 
case the positional coordinates of the spatial point are 
fitting parameters. Preferably, LOR's whose distance of 
closest approach is greater than a preselected distance 
are discarded. Alternatively, the figure of merit is 
generated by obtaining the crossing point or the distance 
of closest approach of each pair of LOR's and calculating 
the standard deviation of the crossing points or the 
obtained distances. 

In accordance with yet another aspect of the 
present invention, a method of PET imaging is disclosed. 
Coincidence radiation events from a calibration source are 
detected with at least two detector heads. Correction 
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factors that correct for mechanical misalignment of the 
detector heads are calculated from the coincidence 
detected calibration source radiation. During a 

diagnostic imaging procedure performed on a subject, image 
data are generated in response to radiation collected with 
the detector heads. The image data are corrected with the 
correction factors. The corrected image data are 
reconstructed into an image representation. 

In accordance with still yet another aspect of 
the present invention, a coincidence imaging system is 
disclosed. The system includes a gantry. A plurality of 
flat panel detectors are disposed about the gantry. A 
data memory stores measured data about radiation events 
detected by the detectors. A calibration memory stores a 
plurality of calibration parameters for correcting data 
measured during a patient scan. A processor in 

communication with the calibration memory and with the 
data memory calculates the calibration parameters by a 
minimization algorithm that includes optimizing fitting 
parameters with respect to acquired radiation data 
associated with a point radiation source. 

Preferably, the calibration parameters include 
parameters relating to positional coordinates of the 
plurality of detectors. The gantry is preferably 
rotatable. The figure of merit is preferably generated by 
summing the square of the distance of closest approach of 
each LOR to a spatial point, in which case the positional 
coordinates of the spatial point are fitting parameters. 
Alternatively, the figure of merit is generated by 
obtaining the crossing point of each pair of LOR's and 
calculating the variance of the crossing points. To 
reduce noise, the minimization algorithm preferably 
discards measured data about radiation events whose energy 
is outside a preselected energy range. 

One advantage of the present invention is that 
the calibration is performed in coincidence imaging mode 
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and does not require mechanical re-configuration the 
system between calibration and PET imaging. 

Another advantage of the present invention is 
that it simultaneously calibrates all three detector head 
positional coordinates, e.g. the radial, tangential, and 
orientational coordinates. 

Another advantage of the present invention is 
that it provides improved image quality. 

Yet another advantage of the present invention 
is that it is compatible with both gamma-PET and dedicated 
PET systems . 

Still further advantages and benefits of the 
present invention will become apparent to those of 
ordinary skill in the art upon reading and understanding 
the following detailed description. 

Brief Description of the Drawings 

The invention may take form in various- 
components and arrangements of components, and in various 
steps and arrangements of steps. The drawings are only 
for the purposes of illustrating preferred embodiments and 
are not to be construed as limiting the invention. 

FIGURE 1 shows a diagrammatic representation of 
a gamma-PET system; 

FIGURE 2A shows a diagrammatic representation of 
a measurement of a point source using a well-aligned 
three-head PET system; 

FIGURE 2B shows a diagrammatic representation of 
a measurement of a point source using a three-head PET 
system in which one head has a radial misalignment Ar; 

FIGURE 2C shows a diagrammatic representation of 
a measurement of a point source using a three-head PET 
system in which one head has a tangential misalignment At; 

FIGURE 2D shows a diagrammatic representation of 
a measurement of a point source using a three-head PET 
system in which one head has a orientational misalignment 
A9; 
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FIGURE 3 is a flowchart of a preferred 
embodiment of the calibration method; 

FIGURE 4 shows two radiation detection events 
corresponding to a single electron-positron annihilation 
5 event, along with the spatial coordinate systems used in 
tabulating and analyzing detection data; 

FIGURE 5 is a diagrammatic illustration of the 
calculation of a preferred figure of merit; and 

FIGURE 6 is a diagrammatic illustration of the 
10 calculation of another preferred figure of merit. 

Detailed Description of the Preferred Embodiments 

_ The invention will be described with reference 

O 

*0 to a three-head gamma-PET system. However, it is to be 

^ appreciated that the invention is not so limited, and may 

^ .... 

m 15 instead be applied to other coincidence imaging systems 

|S such as dedicated PET systems, systems with two or more 

heads, and other systems that reconstruct images based on 

51 a line of response. 

m With reference to FIGURE 1, a detector system 10 

fU 2 0 includes detector heads 12, in the instant gamma-pet 

system three heads denoted 12 % , 12 2 , and 12 3 . Each head 
{2 includes a detector processor unit 14 lr 14 2 , which 

identifies the coordinates of a radiation detection event 
on the detector face. As shown in FIGURE 1, the detectors 

2 5 are mounted on a gantry 16 at 0*, 90*, and 27 0* 

orientations. However, other numbers and orientations of 
heads can be used. The specified orientations are to be 
taken as approximate values only, because as noted above 
mechanical alignment is not precise. It is also pointed 

3 0 out that the gamma-PET system is a versatile medical 

imaging system which also includes SPECT capability. The 
gantry 16 is preferably rotatable. For nuclear medicine 
imaging, a subject (not shown) is positioned in the 
receiving area located within gantry 16. For calibration, 
35 a point source 18 is positioned in the receiving area, as 
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shown. Of course, for medical diagnostic imaging the 
point source is replaced by an associated patient. 

With continuing reference to FIGURE 1, radiation 
detection events detected by detector system 10 are 
5 collected by a LOR calculating circuit 20. The LOR 
calculator 2 0 includes a coincidence detector 22 that 
determines when two events are within a preselected 
temporal window of being simultaneous. When two events 
are determined to be coincident, the identification of the 
10 detector heads, the coordinates of the radiation detection 
point on each detecting head, and the angular orientation 
of the gantry are supplied to a PET LOR calculator 2 4 that 
calculates a ray or line between the two radiation 
\q detection points. The spatial coordinates of the 

uy 15 radiation detection events in the heads are corrected by 

iTJ correction factors retrieved from a correction memory 26. 

*J3 For example, the retrieved correction factors are offsets 

•j! that a correction processor 28 uses to adjust the apparent 

!GB3 

ii spatial locations of the radiation events in the heads 

y 20 prior to the line of response (LOR) calculation. 

=lj Alternatively, the line of response can be calculated 

'"-4 first and then corrected with the correction factors. 

P In SPECT imaging, the identification of the 

detecting head, the coordinates on the detecting head, the 
25 angular orientation of the gantry and heads, and 
identification of the collimator characteristics of the 
head collimator are communicated to a SPECT line of 
response calculator 3 0 that calculates a trajectory or 
line traversed by the received radiation. Correction 
30 factors from the calibration memory 26 are again retrieved 
and used by a correction processor 28' to correct the 
SPECT trajectory or LOR. 

The acquired and corrected LOR data are 
preferably stored in a data memory or buffer 32. A data 
3 5 reconstruction processor 34 reconstructs an electronic 
image representation from the LOR data stored in data 
memory 32 and stores the resultant image representation in 
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an image memory 36. Portions of the stored image 
representation are retrieved by an image processor 38 and 
converted to an appropriate format for display on a 
monitor 40, such as a video, CCD, active matrix, or other 
monitor. Of course, a color printer or other output 
device may also be used to present the data in a 
convenient format . 

To generate the contents of calibration memory 
26, a point radiation source 18 is disposed in view of the 
detectors 12 v 12 2 , 12 3 . Of course, the point radiation 
source 18 has physical size, which is however preferably 
smaller than the desired calibrated image resolution. LOR 
data is measured for the point source 18. The calibration 
memory 26 contents are not applied during this data 
acquisition so that data without any calibration are 
stored in the data memory 32. Alternatively, if 
calibration data from a previous calibration is presently 
stored in calibration memory 26, this calibration may be 
applied as an initial condition (not shown) . A 
calibration unit 42 calculates new calibration data from 
the measurement of the point source 18, and the new 
calibration data is stored in the calibration memory 2 6 
for use in future imaging as described above. 

The method by which the calibration processor 42 
calculates new calibration data is based upon the 
recognition that the LOR's should all intersect at the 
coordinates of the point source 18. When they fail to 
intersect, the calibration processor mathematically varies 
the coordinates of the detector heads until the separation 
of the rays is minimized at a point corresponding to the 
position of the point source 18. It will be particularly 
noticed that the calibration is done in PET mode, with the 
detector system 10 in the PET configuration identical to 
that used for PET measurements. Thus, PET imaging can 
commence immediately after calibration and replacement of 
the point radiation source 18 by a subject, e.g. a 
patient, with no adjustments being made to the detector 
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system 10. This approach advantageously avoids the 
potential for introducing new detector misalignments 
during the detector system 10 reconfiguring. 

With reference to FIGURES 2A to 2D, typical 
contents for configuration memory 2 6 will be described. 
FIGURE 2A shows an ideal, perfectly aligned three-head 
gamma-PET system measuring an ideal point radiation source 
18. Three representative lines of response (LOR) 50, 52, 
5 4 are shown, corresponding to three positron-electron 
annihilation events in the point radiation source 18. The 
LOR 50 is detected by detectors 12 2 and 12 3 ; the LOR's 52 
and 54 are detected by detectors 12 1 and 12 3 . Of course, 
the three event pairs represented by LOR's 50, 52, and 54 
are temporally separated so that the detection system 10 
is able to distinguish the events. Because LOR's 50, 52, 
and 54 originate from the point radiation source 18, in a 
perfectly aligned system the three LOR's intersect at a 
point corresponding exactly to the spatial location of the 
point radiation source 18. 

With reference to FIGURE 2B, the case where 
detector 12., is misaligned by a radial error Ar is 
presented. In FIGURE 2B, the dashed box 12 n is the actual 
position of detector, while box 12 1 b is the position of the 
detector erroneously reported to the data acquisition unit 
20. Because of the radial misalignment Ar, the LOR's 52, 
54 of FIGURE 2A which involved detection events on 
detector 12 1 are now incorrectly calculated in FIGURE 2B as 
52 b and 54 b respectively. The three LOR's 50, 52 b , and 54 b 
do not intersect at a single point. If the detector 
system with the radial misalignment of FIGURE 2B is used 
for gamma-PET imaging without calibration, the resultant 
image will be degraded due to the misalignment Ar. 

FIGURE 2C shows the effect of a tangential 
misalignment At. Again, ideal LOR's 52, 54 of FIGURE 2A 
are incorrectly calculated as 52 c and 54 c respectively, 
resulting in image blurring. FIGURE 2D shows the effect 
of an orientational misalignment A9, which produces 
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incorrect LOR's 52 d and 54 d which again result in image 
degradation. Of course, it is to be recognized that 
misalignments of two or even three positional coordinates 
may occur simultaneously in a given detector head, and 
that any or all detector heads may be misaligned. 
Preferably, calibration memory 26 includes corrections at 
least for radial misalignment Ar, tangential misalignment 
At, and orientational misalignment A9 of each detector 
head 12 1 , 12 2 , 12 3 . 

With reference to FIGURES 1, 3, and 4 a 
preferred method performed by the calibration processor 42 
for calibrating a coincidence imaging system which 
includes a plurality of radiation detectors is disclosed. 
As shown in FIGURE 1, the point radiation source 18 is 
positioned in the receiving area of the detector system 10 
which is configured in the PET imaging mode. The system 
measures 60 a plurality of radiation events associated 
with the point radiation source 18. The calibration 
method uses a minimization method to optimize parameters 
of interest. More specifically, initial values are 
assigned 62 to the fitting parameters. The fitted 
parameters are preferably: Ar i , i=l to 3, where Ar- is a 
correction for the radial coordinate of the ith detector; 
Atj, j=l to 3, where Atj is a correction for the tangential 
coordinate of the jth detector; and A0 k , k=2 to 3, where 
A9 k is a correction for the orientational coordinate of the 
kth detector. It will be noticed that there is no 
correction for the orientational coordinate of the k=l 
detector. This is because one detector defines the 
reference orientation. With this reduction, there are 8 
fitting parameters corresponding to the detector 
positional coordinate corrections. In the step 62, 
initial values for these positional coordinate correction 
parameters is typically set to zero. Alternatively, if 
the calibration memory 2 6 contains values from a previous 
calibration it may be preferable to use these prior values 
to initialize the corresponding fitting parameters. 
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The iterative minimization loop begins by 
calculating the lines of response (LOR's) 64. With 
reference to FIGURE 4 which shows the detector system 10 
in simplified form, an electron 80 and positron 82 
annihilate at a spatial point 84 and the annihilation 
emits two oppositely directed photons. One photon is 
detected at a spot 86 on the detector 12 1# while the other 
photon is detected at a spot 88 on the detector 12 2 . 
Detection events 86 and 88 are recognized as simultaneous 
events because the time separation between the two 
detection events is less than a preselected coincidence 
time window. Preferably, the detections are recorded in 
a format such as (h lf x 1 , y 1 , e 1# h 2 , x 2 , y 2 , e 2 , p) , where 
hf identifies the head, (x,., y ? ) identify the detection spot 
on the face of detector h. using the (x, y) coordinate 
system 90, e 5 is the photon energy, and p is the gantry 
rotation angle 92. The LOR is defined with respect to a 
detector- independent coordinate system, preferably with 
respect to the gantry coordinate system (u, v, w) . The 
position of each head in the gantry coordinate system is 
calculable from the radial coordinate r, tangential 
coordinate t, and the orientational coordinate 0. In 
FIGURE 4, the head 12, has 9=90* and head 12 2 has 0=0*. 
The LOR is then defined by its endpoints which are the 
detection points 86, 88 given in gantry coordinate system 
(u, v, w) as: 

u j = (x. + t (p)) cos(0 +p) + r h ^(p) sin(0^ + p) 

^ (i) 

w = (x. + t. (p) ) sin{0 +p) - r ( p) cos (9 +p) 

where j=l,2 corresponding to the two detector heads. It 
will be noticed that in equation (1) tangent t and radius 
r are written as functions or gantry rotation p, thereby 
accounting for a non-circular orbit. For a circular 
orbit, tangent t and radius r are independent of gantry 
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rotation p. The LOR is thus defined by the detection 
point coordinates (u-,, v 1 , w n ) and (u 2 , v 2 , w 2 ) . In 
performing the LOR calculations, the correction parameter 
values of the present iteration are used. That is, the 
5 position of each radiation detector is corrected by the 
present iteration correction values for the radial, 
tangential, orientational coordinates of the detector when 
applying equation (1) . 

In the iterative minimization method, a figure 
10 of merit is minimized. The figure of merit is selected to 
provide a measure of how closely the LOR's come to passing 
through a single point in space corresponding to the point 
radiation source 18. Alternatively, the closest approach 
*J3 of the LOR's to one another defines a presumed source 

^ 15 region, and the figure of merit would then preferably 

Id define the size of this source region, which as the 

=3 minimization reduces the figure of merit would effectively 

l z \ 

]p reduce the source region toward a point. Two figures of 

« merit are presently preferred, although other figures of 

U 20 merit are also possible and fall within the scope of the 

ny invention. 

.2 With reference to FIGURE 5, a first preferred 

y 

U figure of merit is described. A spatial point 100 is 

defined in space by its gantry coordinates (u 0 , v 0 , w 0 ) , and 

2 5 the figure of merit is calculated by summing the distance 

between the point 100 and the closest approach of each 
LOR. For example, in FIGURE 5 a LOR 102 detected by a 
detection point (u 1# v,, , w t ) on the detector 12 v and by a 
detection point (u 2 , v 2 , w 2 ) on the detector 12 3 , has a 

3 0 closest approach d 102 as indicated on FIGURE 5. 

Alternatively, instead of summing the distances, the sum 
may be over the square of the distance, providing a 
conventional least-squares minimization approach. Given 
the two endpoints of a LOR: (u 1# v^, , w^ , (u 2 , v 2 , w 2 ) and 
35 the coordinates of the spatial point (u 0 , v 0 , w 0 ) , the 
point-to-LOR distance d may be calculated as: 



- 14 - 



PKR 2 0679 / 993 



d 2 = 











2 

+ 










2 




U l 


V 0~ 


V l 


2 


v o 


~ V 1 


W 0 




w o 




u o 




+ 




V 2 


~ V 1 


W 2 


~ W 1 




w 2 


~ w l 


U 2 






U 2" 


U l 


V 2~ 


V l 





(2) 



(U 2 ~U 1 ) 2 + (^-V^) 2 + (W 2 ~W 1 ) 



where the (u, v, w) coordinates are all in the gantry 
coordinate system . 

Whether the figure of merit is summed over d or 
d 2 , the coordinates of the spatial point (u 0 , v Q/ w 0 ) are 
preferably included as fitting parameters in addition to 
the eight detector positional correction parameters 
specified earlier, giving an eleven parameter fit when 
using this figure of merit. However, since the gantry 
coordinates of the approximate point radiation source are 
initially known to within the mechanical alignment 
precision, the initial values for coordinates (u Q/ v 0 , w 0 ) 
are typically close to the final fitted values. 

With reference to FIGURE 6 , a second preferred 
figure of merit characterizes the crossing point (in two 
dimensions) or the distance of closest approach (in three 
dimensions) between each pair of LOR's. In the exemplary 
FIGURE 6, three LOR's 202, 204, 206 cross in two dimensions 
at spatial points with coordinates (u^ v 1f w^ , (u 2 , v 2 , w 2 ) , 
and (u 3 , v 3 , w 3 ) . The figure of merit is then the variance 
or standard deviation of the resultant crossing point 
distribution or distance distribution. This figure of 
merit does not require the additional three fitting 
parameters corresponding to the coordinates of a 
mathematical spatial point. However, because the relation 
between every pair of LOR's is calculated, this figure of 
merit is more computationally involved compared with the 
first figure of merit. 

With the fitting parameters defined and 
initialized and a figure of merit selected, any standard 
minimization technique may be employed to adjust the 
fitting parameters 68 and conditionally iterate until 
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optimization is achieved 70 • The preferred optimization 
method is Powell's method, but other methods such as the 
Marquardt least squares optimization algorithm may also be 
employed. When a minimum is reached, the relevant 
5 correction factors such as radial corrections Ar i , 
tangential corrections Atj, and orientational corrections 
A0 k are extracted 72 and stored in the calibration memory 
26. 

During the calibration process, several 
10 additional error-reduction methods are optionally 
employed. A low activity point radiation source is used 
to minimize the rate of "randoms", which are erroneously 
>Baa identified LOR's resulting from random simultaneous 

detection of two different radiation events. That is, two 
M 15 detection events occur with time separation less than the 

tjj coincidence time window, but the photons dxd not in fact 

jfl originate from a single positron-electron annihilation 

'» event. By using a low activity source, the likelihood of 

ii random simultaneous detection events is reduced. 

^ 2 0 Additionally, using the smallest coincidence time window 

jlj possible further reduces the rate of random LOR 

^ detections. Noise is further reduced by only accepting 

Q 

\2 LOR's which fall within a window around the positron- 

electron annihilation energy. Preferably, a window of 

2 5 approximately 3 0% around the 511 keV positron-electron 

annihilation energy peak is used. 

Finally, if the location of the point source is 
known or can be estimated, LOR's with a closest approach 
greater than a preselected threshold are discarded as 

3 0 being likely randoms or other noise events. A threshold 

distance of 10 mm is preferred. In the case of the first 
figure of merit wherein the distance or distance squared 
between the LOR and a spatial point is calculated, the 
threshold is preferably reduced during minimization 
3 5 process as the estimated location of the point source 
becomes more and more precise. 
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The invention has been described with reference 
to the preferred embodiments/ Obviously, modifications 
and alterations will occur to others upon reading and 
understanding the preceding detailed description. It is 
intended that the invention be construed as including all 
such modifications and alterations insofar as they come 
within the scope of the appended claims or the equivalents 
thereof . 



